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Abstract. We study the position estimation of a mechanical oscillator undergoing 
both detuned parametric amplification and continuous quantum measurement. This 
model, which can be utilised to produce squeezed states, is applied to a general 
optoelectromechanical system. Using a stochastic master equation formalism, we 
derive general formulae for the reduction in position uncertainty of one quadrature of 
motion. The filter for extracting the optimal position estimate from the measurement 
record is derived. We also find that since this scheme does not work far into the 
back-action dominated regime, implementing resolved-sideband cooling improves the 
squeezing only marginally. 



PACS numbers: 42.50.Lc,03.65.Ta,71.36.+c 
1. Introduction 

When continuously monitoring the position of a mechanical oscillator at finite 
temperature, three regimes inevitably emerge pp. In the regime in which even the 
oscillator's Brownian motion cannot be resolved, termed here the bad measurement 
regime, the uncertainty of the position estimate is dominated by thermal noise. In the 
classical measurement regime, the thermal Brownian motion can be resolved but the 
zero-point motion cannot, and the uncertainty is then dominated by shot-noise. Finally, 
a measurement strong enough to resolve the zero-point motion results in the back-action 
dominated regime. The border between the second and third regimes, when considering 
force sensitivity, is usually termed the "standard quantum limit" (SQL) [2]. It is here 
that the uncertainty of the position estimate begins to saturate at the ground state 
variance, provided the measurement is efficient and suitably filtered[l]. 

Optomechanical experiments close to and beyond the SQL are presently becoming 
feasible[3j H], enabling feedback[lj or sideband cooling[5j to near the ground state. 
If a spectroscopic or back-action evading measurement scheme is implemented, this 
regime also allows one to perform quantum tomography [6 J or to produce squeezed 
states[7]. Much attention has been paid to such methods which use, at most, modulated 
measurement signals and linear forces. 
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Once near the ground state, nonlinear classical forces present a natural alternative 
way to prepare quantum states of an oscillator. For example, parametric amplification 
can be used to produce quantum squeezing of one quadrature of motion without 
requiring back-action evasion [8]. However, when using ground-state resolving 
measurement in addition to these nonlinear external forces, one can expect the oscillator 
dynamics to be dominated by back-action, thus limiting the effect of the external control. 
As a result, and as we will discuss in detail, a compromise must be made between strong 
measurement (better cooling) and weak measurement (better control). 

Using the model of an optoelectromechanical system, which exhibits the required 
quantum measurement and classical control [HI [10], we examine the effect of parametric 
amplification on quantum measurement (and vice versa) in various scenarios. 
Parametric amplification has already been demonstrated classically in MEMS and 
NEMS oscillators [HJ [121 HB], an d in the quantum regime for trapped atoms [2] and 
has long been considered as an alternative to back-action evasion [T2]. This is due to 
its nearly limitless ability to "pre-amplify" in-phase position fluctuations so that strong 
measurement is unnecessary for ground-state resolution. 

A parametric drive will also squeeze the unconditional thermal noise in the out- 
of-phase quadrature. In the stable below threshold regime, this squeezing does not 
exceed a factor of two. This is also true if the drive is detuned from resonance |1 5 [ IT6] . 
However, detuning correlates the in-phase and out-of-phase quadratures and therefore 
becomes more interesting in the context of position estimation. Measurement of the 
well-resolved amplified quadrature in this case improves the estimation of the squeezed 
quadrature. This idea has been outlined for the first time recently [17], where it was 
shown that using this approach the potential squeezing is limited only by the oscillator Q. 
Critically, this paper focused on fixed sub-optimal quadratures and an analytic solution 
for the squeezing was only derived for the classical measurement regime. Furthermore, 
experimental implementation of this idea would require the filter — also absent from the 
previous work — which extracts the optimal position estimate from the measurement 
record. Here, we provide a complete analysis of the system, with additional focus on 
phenomena related to optomechanics. 

In this paper, we provide general analytic solutions for the conditional variance of a 
parametrically driven oscillator undergoing continuous measurement. We also determine 
the exact form of the filter which produces an optimal position estimate, where the 
uncertainty of this estimate is given by the conditional variance. Since the quantum 
squeezing is ultimately limited by the temperature of the oscillator, it is natural to 
ask whether combining this scheme with resolved sideband cooling allows an increased 
amount of squeezing. To that end we introduce a master equation that includes resolved 
sideband cooling and discuss the potential of this system for producing squeezed states 
in the good cavity limit. 

This paper is organised as follows. In section [21 we describe the model of the system 
and derive equations governing for the expectation values and variances. In section [31 
we examine the unconditional dynamics and steady-state. In section [U we derive the 
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conditional variances and examine the prospects for quantum squeezing. In section El 
we provide the optimal filter for conditioning and in section [61 we examine the effect of 
resolved sideband cooling on the system. 



2. Model 



2.1. Parametric Drive 

Consider a mechanical mode of frequency u m in which the spring constant can be 
modulated at a frequency near 2uj m . For dielectric oscillators, this is usually achieved 
using a high gradient electric field [121 US]; but other methods exist [HI [19]. In position 
and momentum co-ordinates, the interaction Hamiltonian is 

H = ^- + ^-[k + k r sm(oj d t + 29)}, (1) 

where the parametric drive frequency co d = 2{u m + A). 

For a high-Q oscillator, we can instead use as co-ordinates the quadratures X and 
Y, the measurements of which are the outputs Ix and ly of a lock-in amplifier with a 
continuous position measurement as input (see Figure UK). These operators come from 
a rotating wave transformation at frequency uj d /2 



— —x = Xsin(u d t/2) + Ycos(uj d t/2) , (2) 

where, in terms of new creation and annihilation operators 

X = (a + a ] )/V2 
Y = -i(a-^)/V2, 



so that [X, Y] = i and the ground-state variance for each quadrature is 1/2. The 
Hamiltonian can then be written in the rotating wave approximation as[T7] 

H = HAtfa + ih^ (e 2i V - e- 2W cJ 2 ) , (3) 

where x — w m fc r /2fc , is approximately the peak-to-peak frequency modulation. Since 
we chose a rotating frame with respect to the parametric drive, this Hamiltonian looks 
like a resonance at A with a stationary squeezing operator. Note that the phase 9 has 
no effect on the system dynamics but defines the squeezing axes with respect to the 
chosen quadratures. 



2.2. Measurement 

A stochastic master equation is used to model the measurement and damping[20j. 
Initially, we assume a cavity optomechanical system in the bad cavity limit (k ^> u m , 
where k is the cavity linewidth). We can define a measurement rate [i in this limit, 
given by 
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where g is the optomechanical coupling rate, n is the mean photon number and 
x zp f = ^Jh/muj m is the RMS position due to zero-point motion. This parameter ji 
can be interpreted as a coupling rate to a zero-temperature measurement bath, which 
can be compared with the coupling rate to the thermal bath 7. The ratio /i/7 as well 
as the temperature will be used to define the various measurement regimes, as will be 
discussed in section HI 

This kind of position measurement can be decomposed into quadratures. We 
will limit our analysis to the regime fi <C u m , in which case the stochastic master 
equation can be similarly decomposed. The resulting master equation then resembles 
the well-studied model of heterodyne detection used in optical and microwave systems[TJ. 
Additionally, we assume the measurement signal has no thermal fluctuations. An 
observer's expected knowledge of the observable A then evolves as 

d(i) = -U[A,&\)dt + [2^N + /j](V[a^]A)dt+ [2j(N + 1) + fj]{T>[a]A) dt (5) 
+^jl(H[X}A) dWx + ^(H[Y]A) dW 2 . 

where N is the mean bath phonon number, 7 = u m /Q is the intrinsic damping rate, rj 
is the quantum efficiency and dWi and dW 2 are uncorrelated Wiener processes defining 
the residual noise given the measurement results Ix = dQx/dt and Iy = dQyjdt 



dW x = dQ x - y/^i(X)dt (6) 

dW 2 = dQ Y - v^(F)dt . (7) 

The superoperator 

V[d]A = o)Aa - i(a f al + Ad) a) , (8) 
describes the thermal diffusion and back-action, while 

U[a)A = aA + Aa) - (d + a)) (A) . (9) 



describes the noise being introduced to the measurement record. 
2.3. Evolution of Observables 

Applying ([5]) to the quadratures X and Y gives the evolution of the expectation values 

d(X) = [-( 7 +xcos2^)(A > ) - (A-xsin20)(y)]dt + ^j4^V x dW l + ^/4rjjiC dW 2 (10) 

d(Y) = [-(7-xcos20)(Y) + (A + xsin2#)(X)]dt + ^f^iC dW x + ^miVy dW 2 (11) 

where Vx and Vy are the variances in the X and Y quadratures respectively, and we 
define the covariance C = (XY + YX)/2 - (X)(Y). Setting A = \i = 9 = 0, this 
describes an additional damping in X and reduced damping in Y, both proportional to 
the drive strength %■ I n general, the system is stable if x is below a threshold value 
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Figure 1. a) Schematic of a optomechanical system with an electromechanical 
parametric drive, b) Illustration of an unconditional trajectory in X-Y phase space 
with a detuned drive applied near threshold, in the case A 2 > x 2 . The shaded ellipse 
denotes the conditional uncertainty at a given time based on previous measurements. 
The angles 9 — 7r/4, ao and ai are defined between the Y axis and respective dotted 
lines. 



The detuning causes a continuous rotation between amplified and damped quadratures, 
thus enabling a stronger parametric drive without self-oscillation. For x 2 < A 2 , the 
oscillator's trajectories through X-Y phase space are elliptical with frequency 

(13) 



0J„ 



X 2 , 



and for xlh > X 2 > A 2 , the trajectories take the form [ci smh(iu e ) + c 2 cosh(ic<; e )]e~ 7t . 

We can also find the evolution of the variances by applying the master equation to 
X 2 , Y 2 and XY + YX, then using 

dV A = d(A 2 ) - 2(A)d(A) - (d(i» 2 

dC = -d(XY+YX)-(X)d(Y)-(Y)d(X)-d(X)d(Y) . 

Applying the Ito rules (i.e. (alt) 2 = dWolt = and (alW) 2 = alt) and setting third 
and higher order cumulants to zero (i.e. Gaussian input states) produces quadratic 
differential equations for the variances 

jV x = - 2( 1 + xcos(26))V x - 2(A- X sin(20))C + 1 (2N+1)+^ - 4 W (V 2 + C 2 ) (14) 
jV Y = -2( 7 -xcos(2^))Ky + 2(A + xsin(2^))C + 7(2^+1)+^ -4^(V^ + C 2 ) (15) 
jC = -2 1 C-A(V Y -V x ) + X sH^)(Vx + V Y )-4 w C(V x + V Y ). (16) 



3. Unconditional Steady-State 



Let us for now take the limit where the efficiency rj = 0, implying that all of the 
measurement results are discarded. This gives us the unconditional dynamics, and the 
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variances obtained are those that would be inferred from a spectrum analysis over an 
infinite time. In this case, the measurement strength only appears as an additional 
phonon number due to back-action 

Nba = ^- (17) 

Applying detuning shifts the angle of the squeezing axes significantly. In anticipation of 
this we set 6 = 7r/4, which without detuning will amplify fluctuations maximally along 
an axis rotated through an angle 7r/4 from the Y quadrature (Note: this will henceforth 
be called the antisqueezing angle). Solving ( jT4lfT6i) with the time derivatives equal to 
zero results in the steady-state variances 

v x = (i--r^y (V T + N BA ) (18) 

Vy=(l+ /j^*^ ) (Vr + N BA ) (19) 

C = — (V T + N BA ), (20) 

where Vr = N + 1/2. Clearly, for A « x, the antisqueezing axis is close to Y as it would 
be for 9 = A = 0. 

The exact squeezing is characterized by minimizing the variance over an angle a 
where 

V a = Vx c °s 2 ot + Vy sin 2 a — 2C cos a sin a . (21) 

As long as Vy > Vx, the maximally squeezed and antisqueezed quadratures (VI and 
V + ) become 

V± = ~[(Vx + V Y ) ± (V Y - V x ) sec(2« )] where (22) 
a = - tan' 1 ( —^M—-] , (23) 



2 \VY(t)-V x {t) i 
where a is the antisqueezing angle (see Figure [lb). Applying this to Eqs (I18H20I) gives 
y _ = V T + Nba 
1 + X/Xth 

V + ^ (25) 
1 - XlXth 

a =^tan- 1 (^-), (26) 



where the threshold value Xth, given by ( 1121) defines the maximum drive strength before 
the system self-oscillates. Therefore, for both detuned and resonant drives, a maximum 
unconditional squeezing VL/Vr of — 3dB can be achieved. Squeezing that surpasses this 
amount has been reported in the literature [121 HH], however this was measured using 
only frequency components of the motion near the mechanical resonance. When only 
the resonance is included, the maximum squeezing can be confirmed by using Langevin 
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equations in the Fourier domain to be 6dB[2I]. This is not in contradiction with our 
time-domain approach, which is effectively an integral over all frequency components. 

Note that the squeezing we have discussed in this section is relative to the thermal 
variance, Vl/Vy, which will for our purposes be called classical squeezing. To achieve 
quantum squeezing, the variance must be below the Heisenberg limit. The quantum 
squeezing ratio is then V-/V g , where V g in our units equals 1/2. For quantum squeezing 
of the unconditional variance, the following condition is required 

N + N BA < 1/2 , (27) 

which places a strict upper bound on the measurement strength and temperature. 



4. Conditional Steady-State 

We will now demonstrate that if rj > 0, using a detuned parametric drive enables a 
much greater degree of classical and quantum squeezing. This is because the squeezed 
quadrature can be inferred from previous measurements of the amplified quadrature. 
The conditional variance arising from the master equation can be thought of as the 
mean-square difference between the quadrature amplitude X the observer's optimal 
estimate X est 

V x = ((X(t)-X est (t)) 2 } , (28) 

and similarly for Y. This assuming that X est and Y est are calculated in the correct way 
from the measurement record. We will derive this optimal filter in the next section. 

With no parametric drive, the steady-state conditional variance in both quadratures 
can be obtained by simply setting dV/dt = in ( ITHfTBJ) 



where 



Vo = ^V^- 1 , (29) 

2 

SNR = 2 W {V T + N BA )h = 2 VKN + l/2)/ 7 + ^- , (30) 

is a ratio of mechanical signal to shot-noise. In the classical regime of large N, this can 
be thought of as the ratio of mean-square thermal displacement (Sxt) 2 to (5x 7 ) 2 the 
square of the minimum distance resolvable over a time At ~ 1/ (4rf) pQ 

Taking SNR <C 1 in Eq. (|29|) produces the bad measurement regime, at the limit 
of which Vq reduces to the unconditional variance Vr- As SNR increases past unity so 
that l/(2r}Vr) < fi/j < 2Vt, the conditional variance reduces towards the ground state. 
When /i/7> 2Vr, the back-action in SNR (proportional to (yu/7) 2 ) dominates and the 
strong measurement limit is approached (Vq — > l/(2y/r})). Note that when 2Vr = f] = 1, 
the second regime (i.e. classical measurement) disappears entirely. 

When performing position estimation, a detuned parametric drive results in an 
elliptical gaussian uncertainty. However, the antisqueezing angle of this distribution 
does not in general correspond to a , which defines the axis of the average elliptical 
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Figure 2. (Above) Antisqueezing angle a,\ and (below) squeezing ratio Vx /Vq vs SNR. 
Light curves are for a drive strength \ — IO7 while dark curves are for \ = IOO7. Solid 
lines are for the pump detuning on threshold, while dashed lines indicate detuning 
away from threshold by 7. Note that the squeezing disappears altogether in the strong 
measurement limit. 



trajectory. Figure [lb illustrates this difference, where a parametric drive phase of tt/A 
results in an antisqueezing angle oti for the conditional variance. 

In order to solve the variances for the maximally squeezed and antisqueezed 
quadratures, it is convenient to define the pump phase 9 as a function of SNR and 
other pump parameters so that the X quadrature is always maximally squeezed and the 
covariance vanishes. In order to work in terms of the equivalent antisqueezing angle, we 
make the replacement 9 = 7r/4 — a±. Applying the steady-state condition and C = to 
equations (I14H16P then produces 



v/( 7 + X sin(2a 1 ))2 + 4 7 2 SNR - 7 - x sin(2a 1 ) 

Vx ~ A W (31) 

y/( 7 - X sin(2a 1 )) 2 + 4 7 2 SNR - 7 + X sin(2a 1 ) 
V Y = -. (o/!J 

A(V Y - V x ) 

COS 2«! = 777 —T-r \ ■ 33 

X{V Y + Vx) 

The conditional variances may be solved exactly by finding the antisqueezing angle 



cti in terms of system parameters. The antisqueezing angle (derived in Appendix A ) 
satisfies 




2* j (M) 

From examining the above, we see that the antisqueezing angle increases from the 
unconditional angle «o at SNR = up to 7r/4 in the limit SNR — > 00. Inserting this 
into ( 13TH3"2"|) gives the squeezed and antisqueezed steady-state variances. Notably, since 
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Figure 3. Conditional variance V normalised to thermal variance Vt, for a quadrature 
with no parametric drive (Vo, dotted lines) and in the squeezed quadrature with the 
drive x = IOO7 and detuning at threshold (Vx, solid lines). Initial phonon numbers 
are N — 100 (light), N = 1 (medium) and N = (dark). The increase in squeezed 
variance at high SNR and low temperature demonstrates the effect of back-action on 
quantum squeezing. 



«i does not depend on N, 77 and /1 separately but on their combined form SNR, the 
ratio of squeezed conditional variance to the undriven conditional variance Vx/Vo has 
the same property. These dependences on SNR are plotted in Figure [2j This conditional 
squeezing ratio starts at the unconditional squeezing (rs 1/2), drops to a minimum near 
SNR pa 1, and approaches 1 in the strong measurement limit. 

The degradation of the conditional squeezing for strong measurement occurs even 
in the classical regime where backaction heating is negligible. Once the measurement 
sufficiently strong to directly resolve the squeezed quadrature, the additional benefit 
provided by parametric squeezing is reduced. This can be seen by the squeezed 
and antisqueezed variances (I31H32P becoming independent of the parametric drive as 
SNR — > 00. The effect of the parametric drive in the squeezed quadrature can therefore 
be interpreted as a signal boost from the amplified quadrature, which is of greatest 
benefit near SNR pa 1 where the conditional variance Vq is most sensitive to SNR. This 
is supported by the fact that the squeezed conditional variance Vx starts to significantly 
reduce at a lower SNR than the bare conditional variance V , as illustrated in Figure |3j 

At threshold, Eq. (IMjl reduces to 

/ 2 2 

cos(2a 1 ) = x2 ~ 7 (Vx 2 + 7 2 SNR - 7 v / SNR) . (35) 
In the limit \ 3> 7 and SNR = 1, the following approximation can then be made 

X sin(2«i) pa 7a/2x/7 • (36) 

The amount of squeezing (quantum or classical) achievable is similarly proportional to 
\Jxlli which is limited only by the rotating wave approximation to be less than the 
square root of the oscillator Q factor (i.e. x ^ ^m)- 

At an initial ground state, the optimal regime SNR pa 1 is only on the cusp of 
the back-action dominated regime (assuming the efficiency is near unity). To illustrate 
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Figure 4. Colour plot of the ratio of squeezed variance to ground-state variance for 
X = A = IOO7. A ratio less than unity indicates quantum squeezing. The dashed line 
indicates n/j = N + 1/2, beyond which back-action is dominant. 



the effect of the parametric drive and estimation at low temperature, Vx is plotted 
for various \i and low values of iV in Figure HI Squeezing of the conditional variance 
below the zero-point motion is achieved when Vx < V g = 0.5, which is possible even 
from relatively high initial temperatures or with inefficient detection. The appearance 
of an optimum measurement strength /i at low temperature is in stark contrast with the 
best possible squeezing using a resonant drive, which degrades steadily from a minimum 
0.5V g at /i = to 0.73V g at /z = 7. Even more notably, this scheme vastly outperforms 
back-action evasion in this parameter regime [T7] . 

5. Optimal Filter 

The idea of optimal estimation is based on the fact that since the mechanical noise is 
weighted by a Lorentzian susceptibility, it can be partially filtered out from the uniform 
shot-noise. In real-time control situations such as feedback cooling, the stochastic 
master equation (j3J) can be used to derive a general time- dependent Kalman filter 
for the position estimate [22]. We are interested in a simple filter which produces 
the minimum conditional variances as obtained in the previous section. This can be 
achieved by applying the steady-state condition to the master equation. In this case the 
measurement signals as defined in Eqs (EE) are related to the estimate as follows [23] 

dQx = v / 4^I(X est )dt + dWi (37) 
dQ Y = ^l{Y est )dt + dW 2 . (38) 

Using the above, we can rewrite the mean value equations Eqs (jlOHlip in terms of 
the measurement rather than the white noise residual dW. This results in additional 
damping terms proportional to the measurement rate r//x. As in the previous section, 
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we will choose a pump phase 9 = 71/ A — ot\ and set C = 
Fourier transforming and solving for (X est (u)) gives 



0. Substituting these values 
- xcos(2a ))VVQy(w) 



(39) 

A 2 - x 2 + IxlY + xsin(2ai)(7y - 7*) - lo 2 + iu{j x + 1y) 
where the intrinsic and measurement-induced damping rates of the X and Y quadratures 
are respectively combined as 

j x = 7 + AnnVx (40) 
1y = 7 + Arj/jVy ■ (41) 

This specifies an increasing filter bandwidth as the mechanical signal overtakes the shot- 
noise. These quantities are equal in the bad measurement limit (where rjfi — > 0) as well 
as the strong measurement limit (where Vx ~ Vy ~ 1/(2*/^)), but vastly different near 
SNR w 1 if driven near threshold. We will limit ourselves to the regime A 2 > % 2 , where 
the general form of the above solution in the time domain is 



(X est (t)) = g X xQx(t) * [cos(nt - 
and similarly for the Y estimate 
(Y est (t)) = g YY QY(t) * [cos(nt + 



)e- Tt ] +9xYQ Y (t)* [sm(nt)e- rt ] 



)e~ rt } +g Y xQx(t)*[sm(Qt)e- rt ] . 



(42) 



(43) 



The parameters are obtained from fl39l) and simplified in terms of a% using relations 



from Appendix A 



n 
r 



A 2_ x 2 CO s 2 (2a 1 )(l + 7 2 /A 2 ) 

^(7x + 7y) 
Atan(2«i) 

-1 ( X7cos(2ai) 



tan 



9xx 

9XY 

9yy 
9yx 



sec 4>\J 'ArjfiVx 
A — x cos(2ai) 

n 

sec 0^/4?y/iVy 
A + x cos(2o;i) 



Y 



x ■ 



(44) 

(45) 
(46) 

(47) 
(48) 
(49) 
(50) 



These expressions simplify in the SNR = limit (where cos^ax) = A/xth) to expected 
values (e.g. Q = u e and T = 7). In the high SNR limit (where cos(2ai) = 0) the optimal 
filter has infinite bandwidth and is independent of the parametric drive. 

The measurement records for X and Y, upon applying these Lorentzian filters, 
form an optimal estimate of the oscillator's current position in phase space. In other 
words, the amplitudes of the in-phase and out-of-phase fluctuations are known to within 
uncertainties defined by Vx and Vy. By mixing the estimate back up to a>d/2 and with 
the correct phase shift, an appropriate feedback cooling signal is obtained. In this way, 
conditional squeezing is turned into real squeezing [24J. 
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6. Resolved Sideband Cooling 

Up to this point we have only considered optomechanical systems as ultra-sensitive 
transducers of position fluctuations. It is natural to ask whether there is any 
advantage in using a detuned parametric drive in conjunction with the near- ubiquitous 
optomechanical technique of sideband cooling. The resonant driving case of this has 
been analysed previously, albeit with a focus on the cavity output spectrum[25j. 

We can model sideband cooling by extending the master equation to include a 
cavity bath. This is done by adiabatically eliminating the cavity[26j, and results in the 
back-action terms in equation (JHJ) being replaced by terms analogous to those for the 
thermal bath. The deterministic part of the master equation is then 

(d(i)) = - ^([A,H})dt + 2 1 N(V[a^]A)dt + 2 1 (N + l)(V[a]A)dt (51) 
+ 2 lc N c (V[a)]A) dt + 2 lc (N c + l)(V[a]A) dt , 

where the general forms of the optical damping 7 c and effective cavity temperature Nc 
are given in [5] as T opt and n^. Taking the limit of large cavity loss k 3> u m and zero 
cavity detuning, the back-action noise 7 c(2iVc + 1) from the master equation is equal 
to /i as we have defined in Eq. (jJJ) so our approach is consistent with [5]. 

On the red sideband (i.e. cavity detuning equal to — u m ) and in the good cavity 
limit, the cavity temperature Nc approaches zero and all photons detected are a 
product of phonon absorption. Consequently, only downgoing transitions appear in 
the measurement terms of the master equation. In this regime, adiabatic elimination 
can be performed on the stochastic master equation for the coupled cavity-oscillator 
system as in [22]. The resulting stochastic master equation for the oscillator alone is 
then equivalent to heterodyne detection of a cavity output [23] and has the form 

d(A)= -U[A,H})dt + [2jN}(V[tf}A)dt+[2j(N+l)+v}(V[a}A)dt (52) 
n 

+ vV/2(H[a]A) dW l + ^ W /2(U[ia}A) dW 2 . 

In this case, the Wiener processes are 

dWk = dQ x - y/w(X)dt (53) 
dW 2 = dQ Y - y/w(Y)dt . (54) 

This master equation leads to variance equations by the same method as section |2j 
Letting 9 = tt/A 

jV x = - (2 7 + fi)(y x - 1/2) - 2(A - X )C + 2iV 7 - W [(V X - 1/2) 2 + C 2 ] (55) 
£v Y = - (2 7 + fi)(V Y - 1/2) + 2(A + X )C + 2iV 7 - W [(V Y - 1/2) 2 + C 2 } (56) 

jC = - (2 7 + fiC - A(V Y - V x ) + X(V X + V Y ) - W C(V X + V Y -1). (57) 

Compared to the standard continuous measurement derived earlier, the additional terms 
proportional to /x here are an unconditional linear damping as well as an offset in the 



Position Estimation of a Parametrically Driven Optomechanical System 



13 



1000 



100 



+ 



10 



no cooling 




Ed 





10 



-5 



10 

|x/y 




5 -5 

10 10 



Figure 5. Squeezed conditional variance for an inefficient detector (77 = 0.1), 
normalised to the ground state variance, in the (left) normal and (right) ideal resolved 
sideband regime. In both cases, \ = IOO7 and detuning from threshold is 7. 



conditioning term, such that all measurement terms disappear for a symmetric pure 
state Vx = Vy = 1/2 with C = 0. The threshold condition is now dependent on // 



X f h SB = v% + /i/2) 2 + 



and the unconditional variance with no parametric drive is 
V T 



2 7 ,r 1 

—N + - 

2 7 + /i 2 



(5f 



(59) 



as in [5]. For perfect efficiency, the linear optical damping terms vanish since 

- fi(V x - 1/2) - fi(V x - 1/2) 2 = -fiV x + M /4 , (60) 

leaving variance equations that are identical to (1 1 4ti 1 6 1) apart from an expected factor of 
4 in measurement strength[27j. However, for 77 1, resolved sideband cooling offers a 
qualitative difference since a pure state is always approached in the strong measurement 
limit. Without cooling, this limiting variance increases by a factor of 1/y/rj. 

These differences are small, however, in the important regime around /i ~ 7. This 
is shown in Figure |5j using numerical solutions to the steady-state variance equations. 
With resolved sideband cooling, the squeezing extends further into the back-action 
dominated region fi > 7, however the maximum squeezing is unchanged. Outside 
the good cavity limit, the cavity bath temperature increases and the result can be 
expected to be worse. Therefore, while resolved sideband cooling aids in resolving zero- 
point fluctuations in the presence of detector inefficiency, the requirement of strong 
measurement precludes it from being significantly useful in the context of parametric 
squeezing. Sideband cooling is, however, quite compatible with backaction evasion-based 
squeezing schemes as they both operate best in the (would-be) back-action dominated 
regime and can both be implemented using similar techniques [1]. 
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7. Conclusion 

Unlike classical systems, the strength of any continuous measurement of a quantum 
system plays an important role in its dynamics, with the optimal strength depending 
on the application. Many uses are being found for weak measurement [28] . while 
quantum non-demolition techniques such as back-action evasion work best in the strong 
measurement limit [7J [29] . We have described one case in which a measurement strength 
near the standard quantum limit is preferred, as it is in gravity-wave detectors for similar 
reasons [2]. 

We have shown that with optimal estimation, a detuned parametric drive greatly 
reduces the uncertainty of one quadrature of motion when the signal-to-noise ratio is 
of order unity. At low temperature, since the back-action in this regime is weak, the 
parametric drive allows uncertainties well below the ground state. To our knowledge, 
this is the only steady-state scheme for parametric squeezing for which the squeezing is 
limited solely by device parameters and not by fundamental constraints. 

We have derived the optimal filter for producing the position estimate, which 
can be used to perform feedback and achieve single quadrature confinement to below 
ground-state uncertainty. Confinement based on resolved-sideband cooling only has a 
significant effect in the back-action resolved regime, which is not conducive to squeezing 
by parametric methods. 
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Appendix A. Derivation of Conditional Squeezing Angle 



With the steady-state and C = conditions, the quadratic equations ( TI41IT5]) can be 
expressed as 

2 W V* + (7 - xsin(2ai))W = 2 W V X + (7 + X M^iWx = l{V T + N BA ) . (A.l) 

Rearranging the first equality, 

2 W (V Y + V x )(Vy - V x ) + 7 (W - V x ) = X sin(2 ai )(W + V x ) . (A.2) 

Dividing this through by Vy — V x and using Eq. ( )33|) gives a new equation for the 
antisqueezing angle 

A tan(2ai) = 2rm{V Y + Vy) + 7 , (A.3) 

which as expected, reduces to the unconditional result Eq. (j26j) when 77// = 0. We would 
like a form for a± in terms of system parameters only, for which we can rearrange Eq. 
(lA.ip again in the form 

2 7 (Vr + AM = 2^(1^ + V x ) + 7 (Vy + Vx) - Xsin(2« 1 )(W - V x ) ■ (A.4) 

Dividing (1A.2I) through by Vy + Vx, then substituting for xsin(2ai) above leaves a 
drive-independent relation between the two quadrature variances 

2V Y V X ( W + 1 ^ ) = 7 (V T + N BA ) . (A.5) 



Vy + Vx 

Substituting the relation (from Eq. ( 133]) ) 

l--cos (2«0= (V r y + V3f)a , (A.6) 

as well as Eq. ( 1A.3I) into ( 1A.5I) . we can obtain an expression containing only a function 
of ai, the thermal variance Vr and other parameters. The general form of ol\ is now the 
solution to the equation 

(A 2 tan 2 (2 ai )- 7 2 ) - ^ cos 2 (2 ai )^ = 8 W (F T + N BA ) , (A.7) 



or written in terms of SNR, 



A 2 tan 2 (2 ai ) - x 2 sin 2 (2ai) + cos 2 (2ai) = 7 2 (1 + 4SNR) . (A.8) 

This can be solved analytically by multiplying through by cos 2 (2a!) and transforming 
sine to cosine, then solving a quadratic equation, resulting in 



A /\4 + X 2 +47 ! SNR- v '( x f k - x 2)2+8(xi + x 2 )7 2 SNR+167 4 SNR 2 

cos2ai= d .!\- 1 
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